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Recent photoemission experiments have shown strong surface modifications in the spectra from 
vanadium oxides as (V, Cr)2C>3 or (Sr, Ca)V03. The effective mass is enhanced at the surface and the 
coherent part of the surface spectrum is narrowed as compared to the bulk. The quasiparticle weight 
is more sensitive at the surface than in the bulk against bandwidth variations. We investigate these 
effects theoretically considering the single-band Hubbard model for a film geometry. A simplified 
dynamical mean-field scheme is used to calculate the main features of the interacting layer-dependent 
spectral function. It turns out that the experimentally confirmed effects are inherent properties of a 
system of strongly correlated electrons. The reduction of the weight and the variance of the coherent 
part of the surface spectrum can be traced back to the reduced surface coordination number. Surface 
correlation effects can be strongly amplified by changes of the hopping integrals at the surface. 



I. INTRODUCTION 

A large class of transition- metal oxides display a metal- 
insulator transition (MIT) upon variation of pressure, 
temperature or chemical doping. According to effec- 
tive one-particle theories these materials are normal met- 
als. The observed transitions are caused by correlations. 
Typical materials that show such a behavior are vana- 
dium oxides. There are numerous experimental^ and 
theoretical^ studies which describe and try to explain 
the physics of these materials. The electronic structure 
in the vicinity of the transition has been revealed by pho- 
toemission experiments^. They show a certain internal 
structure of the 3d-derived bands on the metallic side of 
the transition. Essentially, there are two peaks in the 
angle-integrated spectral function A(uo): One is located 
at the Fermi energy and is usually called " coherent part" 
of the spectrum. The other structure , called " incoherent 
part", is located at higher binding energies. The metal- 
insulator transition is then driven by a redistribution of 
spectral weight from the coherent to the incoherent peak. 
At the critical point the weight of the coherent peak van- 
ishes and the system undergoes the transition to the in- 
sulating phase. 

While an effective one-particle theory as the local- 
density approximation within density-functional theory 
(DFT-LDA) can reproduce the coherent part of the spec- 
trum apart from a renormalization factor, it fails to de- 
scribe the incoherent part, which is caused by strong elec- 
tron correlations. On the other hand, the correlation- 
induced MIT is a classical subject of many-body theory. 
If there is no long-range magnetic order in the metal- 
lic as well as in the insulating phase, the transition is 
known as the Mott-Hubbard transition^. For band mod- 
els there are early theories that describe the incoherent 
or the coherent part of the spectrum. Nowadays, with 
the dynamical mean- field theory (DMFT)L&2 there is a 
general approach at hand which includes and unifies both 
aspects. The DMFT reproduces the coherent (usually 
called quasiparticle resonance) as well as the incoherent 
(usually called lower Hubbard band, LHB) part of the 
spectrum. There is a third feature, the upper Hubbard 



band (UHB), which is not seen in photoemission experi- 
ments since it is located well above the Fermi energy. The 
DMFT yields a MIT as described abova&ifi, For bulk 
systems LDA+DMFT calculations even can give quanti- 
tative resultsiii*i£ 

Recently, however, a number of studies have shown 
that the photoemission data have to be interpreted 
carefulljiiiii^i^ii^. Especially, the surface sensitivity of 
photoemission has been reconsidered. It has been shown 
that a proper interpretation of the data cannot rely on 
the presumption that the electronic structure at the sur- 
face is almost the same as in the bulk. In fact, by com- 
paring results from measurements with different surface 
sensitivities, it has been shown that the spectral function 
can be strongly modified at the surface compared with 
the bulk. It has been the intention of these studies to ex- 
tract reliable bulk properties from the experimental data 
which can be compared with the results of LDA+DMFT 
theory, for example. However, also the layer dependence 
and surface modifications of the spectral function have 
been measured. These are interesting on its own. A 
number of qualitative trends can be deduced safely: 

• The weight of the coherent part of the spectral func- 
tion is reduced at the surface. 

• The width or more accurately the variance of the 
coherent peak tends to be somewhat smaller at the 
surface. 

• Changes of the bandwidth (and thus of the effective 
correlation ~) affect the weight of the coherent part 
of the spectral function much more at the surface 
than in the bulk. 

It is an interesting theoretical task to explain these trends 
and to reproduce the respective features in the spectral 
function A(u>). 

In this paper we consider a system in a film geome- 
try which exhibits a Mott-Hubbard MIT and investigate 
the layer dependence of the spectral function. The inter- 
esting question is whether the experimental findings are 
inherent properties of a film geometry of correlated elec- 
trons or whether additional modifications at the surface 
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are required, such as surface relaxation, surface phase 
separation or surface reconstruction. If the measured 
trends are generic they should be present in any model 
system that shows a Mott-Hubbard transition and has a 
layer geometry with surface (s). Thus, though the inves- 
tigated materials usually require the study of multi-band 
models to account for the band degeneracy, the surface 
modifications and layer dependence of the spectral func- 
tion can be investigated qualitatively using a single-band 
Hubbard model in a film geometry. The advantage is that 
a (semi-) analytical approach can be used which makes 
the physical mechanisms beyond the surface effects most 
explicit. 

Of course, such a simple (semi-) analytical theory has 
to retain the fundamental physics of the Mott-Hubbard 
transition as described above. This rules out theories as 
static mean-field theory or Hubbard-III approximation 5 
since they do not give the correct three-peak structure 
of the spectral function. The three-peak structure of 
the spectral function or equivalently a two-peak struc- 
ture of the imaginary part of the self energy is crucial 
for a qualitatively correct (mean-field) description of the 
Mott-Hubbard transition. A recent study^ has shown 
that experimental data can nicely be fitted once this 
condition is met. Here we will use the recently devel- 
oped "two-site" DMFTii. This is an approach which 
keeps the essence of the DMFT but simplifies the mean- 
field equations by a physically motivated approximation 
to allow for analytical calculations. The two-site DMFT 
is an extension of the linearized DMFT«i£ which meets 
the minimal condition mentioned above, i.e. leads to a 
three-peak structure of the spectral function with two 
Hubbard bands (incoherent peaks) and a quasiparticlc 
resonance (coherent peak). 

In the next section we specify the model assumptions 
and develop and discuss the two-site DMFT for a film 
geometry. In section III we discuss the reduction of the 
weight and the variance of the coherent spectral function. 
Section IV considers the topic of the enhanced sensitivity 
of the quasiparticle weight at the surface against band- 
width modifications. 



II. MODEL AND THEORY 

We investigate the single-band Hubbard model for a 
film geometry 
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Here a,(3 — 1, . . . , d label the different layers parallel to 
the film surface(s). d is the film thickness. The sub- 
scripts i and j refer to the sites within a layer and run 
from 1 to ./Vii , where 2Vm is the number of sites per layer 
(7V|| — > oo). In the following the on-site energies and the 
interaction strengths are taken to be layer independent 
t°f — to and U a — U. to = defines the energy zero. 



The hopping integrals t"? are assumed to be non-zero 
between nearest neighbors only and is taken to be layer 
independent, too. The nearest-neighbor hopping defines 



the energy unit t 



a/3 



t 



1. Only at the surface the 



intra-layer hopping tjj — tu or the hopping between the 
surface and the subsurface layer t\"j = t\2 may be modi- 
fied {t\\it\% ^ t) to simulate relaxation processes. 

Exploiting the two-dimensional translational symme- 
try, the layer-dependent Green function reads 



G a p, a (k||,w) = lu + fi - <(k||) - E^k^w) 
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kn is a wave vector of the two-dimensional Brillouin zone. 
The dxd matrix £(kn ) is the Fourier-transformed hopping 
matrix 
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Within the DMFT the self-energy is local. The matrix 
£(kii,u;) is diagonal and k|| independent. In the follow- 
ing the spin index a =| , j. is dropped since we are solely 
interested in the paramagnetic phase. The crucial point 
of the DMFT is the mapping of the lattice (Hubbard) 
model onto an appropriate single-impurity (Anderson) 
model (SI AM). The latter is defined in such a way that 
the impurity Green function and the self-energy of the 
impurity model are equal to the on-site Green function 
and the self-energy of the lattice model, respectively 9 . 
The DMFT is exact in the limit of infinite spatial dimen- 
sions D but can be applied as a proper mean-field theory 
to finite-dimensional systems, too. 

The formulation of the DMFT for a film geometry is 
straightforwardi 9 -: The mapping has to be done for each 
layer. Consequently in our case d different impurity mod- 
els fljj^p, one for each layer a, have to be defined. The im- 
purity models can be solved independently for a — 1, . . . d 
but are coupled indirectly by a set of d self-consistency 
relations 19 . As the usual DMFT for a bulk system, the 
DMFT for a film geometry becomes exact for D — > oo. 

Within the two-site DMFTii the mapping procedure 
is strongly simplified. The respective single-impurity 
Anderson models are replaced by models that consist 
of one correlated impurity site and one bath site only 
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This allows for an exact solution of 



imp ' " 2— site" 

the impurity model. There are four parameters in the 
two-site model, two of them are already fixed: the en- 
ergy level of the correlated site ej^ = t a = and the 
interaction strength U a — U. The remaining two param- 
eters, the hybridization between the two orbitals (sites) 
V a and the one-particle energy of the bath site ei"^ , are 
chosen such that the original self-consistency conditions 
are fulfilled in an integral way for the full spectral func- 
tion and especially for its coherent part. A calculation 
completely analogous to the bulk case (see Ref. 17) yields 
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for hybridization V a : 



i/3 



Here the definition 

z a = I 1 - 



dE Q (w) 
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is used. For metals z a is the quasiparticle weight. Its 
inverse z" 1 is the quasiparticle mass-enhancement fac- 
tor. At the critical interaction for the metal-insulator 
transition there is a divergence of the effective mass: 

z a (U c ) = 0. 

Finally, the one-particle energy of the bath site ec is 
obtained from the second self-consistency condition: 

K = n a (6) 

n„ and n a are the particle densities for the correlated 
site in the impurity model H^ p and for a site in the 
a-th layer of the lattice model, respectively. 

Eqs. 10} and © define the parameters of the (two- 
site) impurity models. The latter can easily be solved 
numerically or (for half filling) analytically. The self- 
energy of the two-site model Y^} site {u) is identified with 
the self-energy of the respective layer E a (w). Now at 
once new quasiparticle weights z a are given by Eq. JS}. 
Via Eq. and the spectral theorem one obtains the 
new particle density for each layer n a . Therewith, a new 
set of impurity models can be defined. This circle has to 
be iterated until self-consistency is reached. 

Almost analytic calculations are possible for the sym- 
metric case of the paramagnetic phase at half-filling. Eq. 

is now trivially fulfilled since n d a = n a = 1 due 
to the manifest particle-hole symmetry, which requires 
el = M = if- Furthermore, the two-site problem can be 
solved analytically, which gives the self-energy 
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Thus Eq. |J5J can be evaluated 
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Introducing the coordination number within a layer q and 
between two layers p (bulk coordination number: q+2p), 
the mapping condition can be written more explicitely 

Va 2 = qtaa z a + P*L+1 Z «+1 + Ptaa-l z a-l ( 9 ) 

For a bulk-system, where all quantities are layer inde- 
pendent, the last two equations give a single non-trivial 
solution and one recovers the Brinkmann-Rica& result 
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The critical interaction, however, is different and given 
by the linearized DMFTiS value 



U c = 6t v / 2p + q 



(11) 



For a film geometry, Eqs. © and 10 constitute a 
set of nonlinear algebraic equations which can easily be 
solved by numerical means - even for large systems. Then 
the self-energy E(w) (Eq. Q) and the Green function 
G(k\\,uj) (Eq. J2|) can be calculated and finally also the 
spectral function 
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To make contact with the experimental findings dis- 
cussed above, we have to determine the weight and the 
variance of the coherent part of the spectrum. These 
quantities can be obtained analytically: The coherent 
peak is described by the "coherent Green function" which 
is obtained by expanding the self-energy for small fre- 
quencies 
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and inserting this expression into Eq. J2J. For the self- 
energy JJJ) the parameters a a and b a are given by 
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This yields for the coherent Green function 
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where t is the real-space hopping matrix. The coherent 
spectral function is given by the imaginary part of the 
on-site element of the coherent Green function matrix: 



At h \u) = -hmG ( i Z h J(u + iO + ) 
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Its weight w a and its variance are determined by the 
spectral moments = J° 
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To obtain the spectral moments, we perform a high- 
frequency expansion of the coherent Green function (see 
Ref. 20) 
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FIG. 1: Layer-dependent spectral function A a (u>) of a 21 layer 
sc(100) film for uniform hopping t = 1. 

This yields 

w a = z a 

A 2 a = z a (gtl a z a +pt 2 aa+1 z a+ i+pt 2 aa _ l z a -i)(19) 

With these results at hand we can analyze the surface 
effects of the coherent part of the spectrum. Eqs. (JSJ 
and JHJ turn out to be sufficient to discuss most surface 
effects. Additionally, we can calculate the full spectral 
function using Eq. H12|) . 

III. SURFACE WEIGHT AND VARIANCE 

Recall that the experiments show three characteristic 
surface modifications of the spectral density: a reduced 
weight and a tendency to a reduced variance of the co- 
herent part of the spectrum and an enhanced sensitivity 
against bandwidth modifications at the surface. Let us 
start with the first two phenomena, which have been ob- 
served for (La, Ca)V0 3 (Ref. 13 ), for (Ca, Sr)V0 3 (Ref. 
14 and Ref. 15) and for (V,Cr) 2 3 (Ref. 11). We will 
discuss the spectra from CaV0 3 and SrV0 3 (Ref. 14, 
15) as representative examples. In the spectra of both 
materials there are two peaks, the coherent peak near 
the Fermi energy and the incoherent peak around 1.5 eV 
below Ep. This general structure is found for both, a 
surface-sensitive measurement as well as for a measure- 
ment with a weak surface sensitivity. The weight of the 
coherent peak is considerably reduced at the surface and 
its variance is somewhat smaller. The question we like to 
address is the following: Is the surface weight and vari- 
ance reduction solely caused by the reduced coordination 
number at the surface or are these effects rather caused 
by modifications of surface parameters (which may be 
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FIG. 2: Coherent part of the layer-dependent spectral func- 
tion from Fig. Q (solid line), the free spectral function (U = 
dashed line) and the free spectral function narrowed by a con- 
stant factor 2 = 0.46 (dotted line). 



due to e.g. surface relaxation or reconstruction). To 
this end we compare results for a Hubbard film with uni- 
form hopping t with those obtained from calculations for 
tii, ti2 7^ t at the film surface. 

Uniform parameters. We start the discussion with 
the case of uniform parameters. Fig. shows the layer- 
dependent spectral function A a (uu) of a 21 layer simple- 
cubic (100) (sc(100)) film for U = 10 and half-filling. 
Layer "1" denotes the surface layer. The central layer 
("11") simulates bulk properties rather well. Each spec- 
trum consists of a coherent peak around the Fermi en- 
ergy and incoherent peaks at u> ~ ±8. In the photoemis- 
sion experiments only the occupied part of the spectrum 
is seen which leads to the observed two-peak structure. 
Compared with the experiments the different peaks are 
separated more clearly in Fig. ^ This is due to the fact 
that damping effects due to a finite imaginary part of 
the self energy are neglected completely in the two-site 
DMFT. Damping effects are expected to be less impor- 
tant for the coherent peak shown in Fig. [21 As can 
be seen in the figure and as is also noticed in Ref. 15, 
the shape of the coherent spectrum is well described by 
the uncorrelated (U — 0) spectral function, apart from 
a correlation-induced scaling factor. This factor z can 
be identified with the quasiparticle weight z from Eq. 
(JSJ). In Fig. [21 comparison is made with the free spectral 
function narrowed by a constant factor z — 0.46 for all 
layers (dotted lines). While this factor works well for the 
surface it is too small for the bulk. Hence the scaling 
factor and therefore the quasiparticle weight is reduced 
at the surface. This is the same trend as seen in the 
experiments. 
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FIG. 3: z as a function of U in the bulk and at different sur- 
faces of a simple cubic lattice for uniform parameters. Inset: 
the ratio of the surface and the bulk variance of the coherent 
peak as a function of U. 




FIG. 4: zi and zn ~ Zb as functions of U for a sc(100) film 
(d — 21) and reduced hopping parameters at the surface. 
Inset: The ratio of surface and bulk variance of the coherent 
peak as a function of U. 



However, the reduction is much weaker as in the exper- 
iments. A possible reason is that the sc(100) surface is 
rather closed, i.e. the ratio between the surface and the 
bulk coordination number %p+ q = f ^ s near unity. In Fig. 
we compare the quasiparticle weight z of the sc(100) 
surface with the ones of more open surfaces (sc(110) and 
sc(lll)). Indeed the surface quasiparticle weight is more 
reduced if the surface is less closed. 

In the inset the ratio between the surface and the bulk 
variance is shown. The surface variance is reduced for all 
interactions U which is in agreement with experiment. 
Again the reduction is more pronounced for the more 
open sc(lll) surface (dot-dashed line). We can distin- 
guish between a direct effect that is already present at 
U = and is given by the ratio between the surface and 
the bulk coordination number (see Eq. (|19|l ) and an ad- 
ditional indirect effect for finite interactions U which is 
due to correlations. 

The results discussed so far are not at all specific to 
the considered low-index surfaces of a simple-cubic model 
structure. On the contrary, one can show analytically 
that (within the two-site DMFT) the reduction of the 
weight and the variance is inherent to any Mott-Hubbard 
system, irrespective of the special geometry. We consider 
a semi-infinite system (d — > oo) with uniform parameters. 
We will show that the weight at the surface Z\ is smaller 
than the weight in the bulk Zi, (0 < Z\ < Zb < 1 for 
< U < U c ). From Eqs. © and © one readily derives 
the recursion: 



for the surface [a = 1, zq = 0) 



2i = l- 
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(21) 



On the other hand, from (|10|1 and we have for the 
bulk 



Zb = 1 



U' 2 
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(22) 



Now suppose for the moment that for a given U (0 < 
U <U C ) 



Z\ > z b 



(23) 



We will show that this leads to a contradiction. It follows 
from (|21|l and l|22|) that the assumption (|23|l requires 



Z 2 > Zl> z b . 

Now assume that the relation 

Z a > Z Q _l > z b 



(24) 
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is valid for arbitrary (fixed) a. Then j_^ b > 1 holds. We 
immediately get 



1- Zb ( 1 - z h z a -i \ q ( 1 - Zb 



1 Z a \ 1 Z a Za J p \ I Z a 



- 1 > 1 (26) 



z a +i — z a — z Q _i, (,zuj 



since the last two addends are positive. After some alge- 
bra one arrives at 



where U' = ¥ is defined for convenience. Now we find 



(1 -z b )(2p+q) q z a _ 
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Using 12211 and multiplying the inequality with z a gives 



z a U' q 
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The left hand side of this expression is just z a+ i as seen 
from Eq. lf2*U|) . Hence 



1 > Za > Zb 



(29) 



It is shown that (|2*9"jl directly follows from lj2*5l) and (|2*5|) 
holds for a — 2 (|2~ljl . Consequently (|2*9"|l holds for all a. 

On the other hand z a has to converge against the bulk 
value zt for a — > oo 



(30) 



The last two statements (|29|) and <|3U[) are contradictory. 
Hence our assumption (|23|l . i.e. that 0i equals or is 
greater than Zb, cannot be valid. Thus it is shown that 
the quasiparticle weight is reduced at the surface 



z x < z b for < U < U c 



(31) 



for uniform parameters, irrespective of the geometry. 

A short and vivid explanation of this finding is based 
on the reduced coordination number at the surface. Let 
us compare Eq. J§J) for the surface and for the bulk. 



Vt = ziq + z 2 p , 
V b 2 = z b q + z b p + z b p 



(32) 



There is a (positive) addend missing in the surface ex- 
pression. This tends to reduce the surface parameter V\ 
which causes a reduced quasiparticle weight z\{U) in turn 
(Eq. ©). 

All experiments discussed above show a reduced weight 
of the coherent peak at the surface. Some of them may 
even be interpreted as showing an insulating surface of 
a metallic bulk. However this is ruled out for uniform 
parameters. From Eq. 12U|) it follows immediately, that 
all layers are insulating if the quasiparticle weight of the 
first layer vanishes. (Recall that zq = 0.) 

Similar considerations as for the weight apply for the 
variance. From Eq. (1 191 and © we have: 



Aj = ziVi 



z b V b 2 



(33) 



Using Vl = Yfr a U ( from E( l- ©) and E( l- © we 
end up with 
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(34) 



Since we already know < z\ < Zb < 1 this immediately 
yields 

Al<A 2 




FIG. 5: As Fig. 0but for reduced hopping parameters at the 
surface. 



for any surface geometry and for uniform but arbitrary 
parameters. 

Again, a closer look at JSJl and 10 gives a vivid expla- 
nation of this finding: Starting from (|19|) . we get for the 
surface and the bulk variance: 



A\ = zi(qzi +pz 2 ) , 
Al = z b (qz b + pz b + pzb) 



(35) 



Here one can see both effects that have been already dis- 
cussed at Fig. |3J The direct effect is just the missing 
addend in the surface expression (similar to l|32|l ). This 
reduces the surface variance even for U = 0. The indi- 
rect effect, which is due to correlations, is caused by the 
reduced surface weight zi(U) < Zb{U). Both effects add 
and result in a reduced variance at the surface as seen in 
Fig. El 

Modified surface parameters. Even for the rela- 
tively open surface sc(lll) the reduction of the quasi- 
particle weight at the surface is smaller as compared to 
the reduction found in the experiments. We therefore 
want to investigate whether a modified hopping at the 
surface (caused e.g. by a surface relaxation) can further 
reduce the spectral weight and the variance of the coher- 
ent peak. Indeed Fig. 0] shows a considerable influence 
of the surface hopping on the quasiparticle weight. If 
both, the hopping within the surface layer tn and the 
hopping between the surface and the subsurface layer t\2 
are reduced (dot-dashed line) the weight of the coherent 
peak is almost negligible for a wide range of interactions 
strengths. This applies e.g. for U = 10 which is still far 
away from the bulk metal-insulator transition. Never- 
theless, looking at the layer dependent spectral function 
(Fig. |3J), there are strong correlation effects at the sur- 
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FIG. 6: The slope of the quasiparticle weight as a function of 
the hopping parameter t for constant U = 10. 



face. Now there is only a very narrow peak at the Fermi 
energy for a — 1. This fits well to the experiments cited 
above. The inset of Fig. 0] shows analogous trends for 
the variance. 

Summing up, one can state that both experimentally 
established trends, i.e. the reduced weight and the re- 
duced variance of the coherent surface peak, are inherent 
features of a Mott-Hubbard system in a film geometry. 
Reduced hopping integrals at the surface may amplify 
both effects considerably. 



IV. SENSITIVITY AGAINST BANDWIDTH 
VARIATIONS 

Another surface effect is mentioned in Ref. 15 for 
(Sr, Ca)VC>3. Again, measurements with strong as well 
as with weak surface sensitivity have been performed. 
It is found that the surface weight z\ is much smaller 
in CaV03 than in SrVC>3 while the quasiparticle weight 
in the bulk zi, is almost the same. Both materials are 
expected to exhibit comparable interactions U. The dif- 
ference is mainly the hopping integral t which is smaller 



in CaV03. In other words, the experiments show that 
the slope of the function z(t) - 4| - is positive and much 
higher at the surface compared to the bulk. Fig. Elshows 
this slope as a function of the hopping t (and thus as a 
function of the bandwidth W = 12t) at a fixed interaction 
U = 10. Below the critical hopping t cr it =0.68 the sys- 
tem is insulating. For the first case, a sc(100) surface, we 
can qualitatively reproduce the experimental result for a 
wide parameter range (t > 0.85). However the enhance- 
ment is very weak. As for the other surface effects, the 
enhancement becomes stronger for more open surfaces 
as well as for reduced hopping integrals at the surface. 
The experimental finding =Q- > ^ is thus expected for 
most parameters. However, the opposite scenario may 
be found as well for smaller hopping integrals t. 



V. SUMMARY 

We have reproduced and explained a number of surface 
effects recently detected by photoemission experiments 
from metallic Vanadium oxides. Within a simplified but 
reasonable DMFT scheme, it was shown analytically and 
numerically that the reduction of the quasiparticle weight 
and the variance of the coherent part of the spectrum at 
the surface are inherent properties of a system of corre- 
lated electrons. It was explicitely shown that these effects 
are caused by the reduced surface coordination number. 
Consequently, the surface modifications are stronger for 
an open surface. These effects can be amplified if the 
hopping integrals at the surface are reduced. 

The situation is not so clear for the enhanced sen- 
sitivity of the surface quasiparticle weight z\ against 
bandwidth variations. Though the experimental trend 



dt 



> 



dt 



is found for a wide range of parameters, the 
opposite is true for parameters that are still realistic. 
Furthermore we did not find any region in the parame- 
ter space where z\ changes considerably but z& is nearly 
insensitive against bandwidth variations. This subject 
needs also more clarification experimentally: In Ref. 14 
for instance, where the same materials are investigated 
as in Ref. 15, the bulk quasiparticle weight still seems to 
change if one compares SrVOz and CaVO^. 
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